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The pressure of hot QCD has recently been determined to the last perturbatively computable order g^\ng 
by Kajantic et al. using three-dimensional effective theories. A similar method is applied here to the 
pressure in the presence of small but non-vanishing quark chemical potentials, and the result is used to 
derive the quark number susceptibilities in the limit /i = 0. The diagonal quark number susceptibility of 
QCD with Uf flavours of massless quarks is evaluated to order g^lng and compared with recent lattice 
simulations. It is observed that the results qualitatively resemble the lattice ones, and that when combined 
with the fully perturbative but yet undetermined g^ term they may well explain the behaviour of the lattice 
data for a wide range of temperatures. 
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1. Introduction 



The grand potential is the most fundamental function describing the equilibrium properties of a thermo- 
dynamic system. By itself it amounts to minus the pressure p{T, fj,i, fj,2, ...) times the volume and from its 
derivatives one immediately obtains such quantities as entropy, specific heats, number densities and suscep- 
tibilities. In quantum chromodynamics a lot of attention has been devoted to determining the pressure in 
the quark-qluon plasma phase and e.g. lattice simulations have been successfully applied to the problem. 
They are, however, only available for temperatures no more than a few times above the critical temperature 
Tc of the deconfinement phase transition and have only recently been extended to finite baryon densities. 

So far the most powerful tool available for analytic calculations in QCD has been perturbation theory, 
the use of which is justified by the small value of the coupling constant g at high energy densities due to 
asymptotic freedom. At = the perturbative series of the pressure has recently been extended to order 
g^lng which marks the final step in a series of impressive computations starting from the determination 
of the order g^ contribution [2] and leading through the orders g^ gHng 0, 5^ (S] and [HIIZI- The 
next 0{g^) term in the series is already out of reach for analytic computations due to infrared problems [Sj. 

At non-zero quark chemical potentials the pressure is at present known only to order g'^lng ^ at finite T 
and to order at T = j9j reflecting in part the computational complications induced by a finite value of /i. 
From these results one may derive the quark number susceptibilities x defined as second derivatives of the 
pressure with respect to the chemical potentials. They are both important and very interesting quantities, 
since they describe the effects of finite density being at the same time directly measurable on a lattice in 
the limit fi (see e.g. ^J^JEl)- The perturbative results for x obtained from ^ in this limit are, 
however, unable to produce even the qualitative behaviour of the lattice data. Recent HTL computations 
(see e.g. ^3]) have improved the situation somewhat. 

In the present paper the pressure of quark-gluon plasma with massless quarks is computed to order g^\ng 
in the limit of small but non- vanishing chemical potentials < /i/T ^ 1. The calculation is a generalization 
of the recent = paper ^ and applies the same method, matching of effective three-dimensional theories 
to four-dimensional QCD , in separating the contributions of different momentum scales to pqcd ■ The 
result is used in deriving the quark number susceptibilities at = to order g^ln g, and the diagonal 
susceptibility is subsequently analyzed and compared with results of lattice simulations. One observes that 
the results of the present computation follow the same trend as the lattice ones, but that the effects of the 
yet undetermined contributions may still affect their behaviour considerably. 

The non-trivial part of deriving the susceptibilities is the evaluation of the fermionic three-loop diagrams 
of full QCD that contribute to the pressure. To do that one applies here the results of yet unpublished work 
|15| . in which the necessary diagrams are computed for arbitrary T and /i. In this paper the results are, 
however, only quoted as expanded to second order with respect to fJ./T, which is sufficient for the present 
purposes. Due to the use of the effective theory approach one may here perform the calculations without 
any resummations applying dimensional regularization in the MS scheme is to regulate both ultraviolet and 
infrared divergences. All fields are considered massless, while the chemical potentials of the quark flavours 
are regarded as being independent and non-zero and the temperature higher than Tc- 

The general setup of the paper is presented in section 2, while the results for the pressure are assembled 
in section 3. Section 4 is then devoted to explaining the computation of the matching coefficients needed 
in deriving the pressure, and the diagonal quark number susceptibility at /i = is discussed in section 5. 
Conclusions are finally drawn in section 6 and the values of the matching coefficients listed in Appendix A. 

2. Setup and notation 

In Euclidean metric quantum chromodynamics is defined by the Lagrangian density 

J^QCD = ^F^,F^,+{pip^lj, (2.1) 



1 



where 

= d^^Al-d^A^+gr^'^AlAl, (2.2) 
= d^~igA^ = d^-igAlT'', (2.3) 

and where the massless quark fields have been combined into a multi-component spinor ip. The Nc x Nc 
-matrices T" are here the generators of the fundamental representation of SU{Nc), and the relevant group 
theory factors read 

(2.4) 
(2.5) 
(2.6) 
and 

DScd = d^'^d"'"' = ^^^^^-d. (2.7) 

The dimensions of the representations are (Ia = 5aa = — 1 for the adjoint one and dp = Su — cLaTf/Cf — 
NcUf for the fermionic one. 

The partition function is defined as a path integral over all fields of the functional 
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where (3 = 1/T and is a diagonal Uf x Uf -matrix in flavour space representing the different chemical 
potentials of the quark flavours. In H2.8|l the space dimensionality is d = 3 — 2e signifying the use of 
dimensional regularization. The partition function gives the pressure through the relation p — T/V\nZ, 
where the infinite- volume limit is assumed. 

At high temperatures and small chemical potentials ji <^T the pressure of QCD can be separated into 
three parts Pqcd = Pe + Pm + Pg corresponding respectively to the different momentum scales 27rT, gT 
and g^T contributing to it ^1[7]. By definition 

Pe(T,a*) = pQCD(T,A*)-^lnyPA^I?Agexp{-5E}, (2.9) 
where is the action of a three-dimensional effective theory with the Lagrangian density 

Ce = iTrf^2 +Tr[A,^o]'+m2TrA2 + ^X^^TrAg + 5/:E, (2.10) 

obtained from QCD by dimensional reduction Tf. In l|2.1()|l niB and the coupling constant gE appearing 
in Fij are parameters to be determined in full QCD, and terms that contribute to the pressure starting at 
order g^ or higher have been assembled in SCe- 

The functions pm (T, fj.) and pQ (T) are similarly defined by 

PM (T, /i) = PQCD (T, ^l)-pEiT,^i)~^lnJ VA1 exp { - 5m } 

= P(^cT,{T,^l)-pE{T,^i)-pG{T), (2.11) 
Cy, = ]^TyFf^+5CM (2.12) 
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with yet another couphng constant appearing in the definition of Fij. At leading order the different 
parts contribute to the pressure as pe ^ 5", Pm ~ 9'^ and pQ ^ g^lng. As is indicated above by writing 
Pq{T), it will be seen that the dependence of pq on fi is of higher order than what is considered here. 

The contributions of the different momentum scales to the /i = pressure have been analyzed in detail 
in T. Since the effects of finite chemical potentials with few exceptions manifest themselves merely as 
changes in the matching coefheients defined there, the treatment of the pressure in the present paper will be 
restricted to quoting the results for the pn's from pP and discussing the effects of finite fj, on the coefficients. 

The momentum integration measure and the shorthands for sum-integrals used from here onwards are 



2tt Y \ in J J {2n 

T yf, (2.14) 



where A is the MS scale and po = 2nnT stands for bosonic and {po} = (2n + l)7rT — ifj, for fermionic 
Matsubara frequencies. In the following sections the chemical potentials usually appear in the dimensionless 
combination 



3. The pressure at < /i < T 



Collecting results from and |16| one may, in analogy with T , write the different parameters required in 
computing pqcD to order g^\n.g in the form 
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where the effects of finite fi can be explicitly seen only in the appearance of an additional term proportional 
to aM2 in (I3.2|l . Here is the renormalized coupling of full QCD, and the values of the matching coefficients 
a can be immediately obtained from the results of 1], |2, |13| and |16) with the exception of ctEs and ckE5- 
While the latter is found almost trivially through a one-loop computation, determining the first one requires 
the evaluation of all three-loop diagrams of full QCD that contribute to the pressure. This procedure is 
explained in section 4. 

Adding together (|3.1|l , H3.2|) and H3.5|) one has obtained a compact expression for the perturbative expan- 
sion of the pressure up to 0{gHiig) 
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with the values of the coefficients a listed in Appendix A and the pole of q;e3 exactly cancelling the 1/e term 
appearing in the order g'^ contribution. One should notice here that the numerical factors appearing inside 
the logarithms of the last g^\ng term can be unambiguously defined only after the full order g^ contribution 
has been determined. 



4. The computation of a^s 

In order to find the matching coefficient Q!E3 one needs to compute the diagrammatic expansion of the QCD 
pressure to three-loop order at finite T and /i, but without any resummations. Since the corresponding 
expansion has already been derived for /i = l^j, one may restrict the treatment here to the /i-specific part, 
i.e. to the diagrams containing fermionic lines. These diagrams are shown in Fig. 3 and contribute to the 
expansion as 

Worm = pLn. + lla + llb + \lc + ^Id + \le + llf + \lg + llh, (4.1) 

where Pj^^^ represents the pressure of non-interacting quarks (see e.g. ^3). 

Applying the usual Euclidean space finite temperature Feynman rules to the individual diagrams, per- 
forming several shifts of integration momenta and taking advantage of the fact that the purely bosonic 
version of 7~ vanishes at 0{e^) 'S', the diagrams may be rewritten in terms of the sum- integrals 

IT ^ fj^^ (4-2) 

in ^ i TtC, (4.3) 
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Figure 1: The two- and three-loop fermionic diagrams contributing to the pressure. The sohd, wiggly and 
dashed hnes stand respectively for the quark, gluon and ghost fields. 
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defined here in analogy with For example, for the diagram b one obtains in the Feynman gauge 
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Substituting Eqs. (|4.8|l - 14.15|l into (|4.1|l and allowing explicitly flavour-dependent chemical potentials. 



one has obtained a representation for pform in terms of the sum-integrals defined above 



Pfc 



Pfcrm 



2 (1 - e) dATpZyi^^ (1° - 22^ 



+ dAg^ [CaTf (1 - e) [8 (1 + e) X^X?!^ - 4e25'T + 425'T 

+ 2e7Wo,o - (2 + e) A^o.o - 2M^_i - 47W_2,2 

- 2CfTf (1 - e) [2 (1 - e) {l° - ' 1° - 4J?T + 4 (1 + e) T^f 
+ 2Mofi - (2 + e) A/'o^o + 2 (1 - e) A4i,_i - 2^.- 

- dAg'nls (1 + e) 2?K] + 4 (1 - e) T^] + 2?^] 



(4.16) 



- 16 [ll[fif]r\pg]+ll[fig]T\fif]) ~ 2(1- e)MoAf^f, fig]- 

The renormahzation coefficient of the gauge couphng, Zg, is given by 

^ 11Ca-4Tf 52 1 
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and in the latter sum of (|4.16() the A/'n,_„'s depend on both /if and /Zg through the respective fermionic 
momenta. The sum-integrals appearing in the result are computed^ in |15| for arbitrary T and /i, and thus 
the coefficient oes is available by simply combining (|4.16l) with the bosonic part of the strict perturbation 
expansion of the pressure (see e.g. Eq. (31) of ITj). The outcome of the computation expanded to second 
order in p. is displayed in Eq. HA.6|) of Appendix A. 



5. Quark number susceptibilities at /i = 

From the expression ()3.8|l one may at once extract the quark number susceptibilities defined by 
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where the symmetry between the massless quark flavours at ^ = has been exploited in the latter equality 
in dividing the result explicitly into a diagonal and an off-diagonal part. The off-diagonal susceptibility x is 
already known to order g^ln g "13", but the result obtained here for the diagonal one x is new. Specializing 
to the physical case A^c = 3 one readily obtains 
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Figure 2: x/xo plotted as a function of T/Tc for n/ = (left) and nf — 2 (right) with various values of C. 
The parameter A is defined as the ratio of the coefficients of the f;^ln g and terms in the expansion H5.2|l . 
The lattice data is from ^UEJj and the values rc/A^|„j,=o — 1-15 and Tc/A^\nf=2 — 0.49 [201 are used. 
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(5.2) 



(5.3) 



has been used to determine the form of the second but last term, which cancels the scale-dependence of the 
lower-order contributions. The symbol p is used in 15.2|l in place of the scale of dimensional regularization 
to emphasize the fact that its value is arbitrary, though one expects it to be of order 2ttT . 

The last term of H5.2|l . proportional to C{nf), represents the yet undetermined 0{g^) contribution to the 
diagonal susceptibility and can only be obtained through a massive computation involving the evaluation 
of all four-loop diagrams of full QCD contributing to the pressure. It is, however, worth noting that unlike 
in the case of the pressure ^ no lattice simulations will be required in this process due to the fact that at 

■^Recently there has appeared a new paper 1191 . in which one- and two-loop sum-integrals analogous to those considered in 
1151 have been independently calculated. 
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Figure 3: x/xo|n/=2 at different perturbative orders for C(2) = —45 or A w —0.18 (left) and the absolute 
values of the individual terms of the series (right). The lattice results are from 

order the contribution of pG to pqcD is /x-independent. It is therefore obvious that C{nf) has no direct 
relation to the coefficient 6 defined in "TI to represent the 0{g^) part of pqcd- 

In Figs. 2 and 3 the result (|5.2I) is plotted as a function of T together with lattice data obtained 
from |10[ Fig. 2 contains plots of the Uf — and Uf — 2 cases for different C{nf) showing that 
for reasonable values of C the perturbative series can be made to approach the lattice results even at 
temperatures surprisingly close to T^. This should, however, by no means be considered an argument 
suggesting that the region of applicability of the present results can be extended down to T^. The success of 
these perturbative predictions at T < 3Tc is in large part merely a consequence of having a free parameter 
available for plotting, and the behaviour of the result at these temperatures may be completely distorted 
by the eventual computation of C(n/). For reference, the susceptibility for n/ = 2 and C(2) = —45 is 
nevertheless plotted in Fig. 3 to different perturbative orders alongside with the absolute values of the 
different terms of the series. For temperatures higher than a few times Tc the convergence properties of the 
series appear to be reasonably good. In the plots the value of the scale parameter p is set to 6.742r for 
Uf = and to 8.112T for n/ = 2, for which the one-loop corrections to vanish |21| . 

A few words about the curious concept of susceptibility at zero quark flavours are probably in order. 
The quantity has been computed above as the formal Uf ^ limit of the general result H5.2|) and has 
been plotted in Fig. 2. The interest in this unphysical limit is due to the existence of recent lattice results 
from quenched QCD |1(J) that provide an interesting and powerful test for the validity of the perturbative 
expansion of the susceptibility. In particular, being able to compare the perturbation theory results and 
the lattice data for different numbers of fermion flavours makes it possible to qualitatively study the Uf- 
dependence of the yet unknown coefficient C(nf). 

6. Conclusions 

In this paper the diagonal quark number susceptibility of QCD at vanishing chemical potentials has been 
computed to order g^\n.g in perturbation theory. This is a three-order improvement to the previous result 
0]. Since the next term in the series requires the determination of the perturbative part of the 0{g^) 
pressure at finite T and and even the corresponding /x = computation seems to be out of reach for 
current computational techniques, the result obtained here will most likely not be subject to improvement 



in near future. Nevertheless, a need for further work aiming at the determination of the 0{g^) terms in the 
expansions of both the pressure and the susceptibilities is obviously present. 

The result obtained here for the diagonal susceptibility resembles the lattice data available but, as demon- 
strated in the previous sections, may also be considerably modified by the yet undetermined contributions to 
X revealing its lack of predictive power. This could actually have been anticipated due to similar problems 
encountered in connection with the perturbative expansion of the pressure. The problem can be viewed 
as a natural consequence of the relatively large value of the QCD coupling constant at the energy scale 
considered here and of the slow convergence of the perturbative expansion of the pressure in the 3d sector. 

While the high temperature, small chemical potential region of the QCD phase diagram is without doubt of 
considerable special interest, it is also worthwhile to study how the results obtained there can be generalized 
to other parts of the (T, /x) -plane. An especially interesting task waiting to be tackled is the building of a 
bridge between the order perturbative results for the pressure at T ^ 0, ^ = ^ and T = 0, fi ^ El- 
In the T fj, region the true ground state of the theory may naturally be modified by the appearance of 
non-perturbative qq-condensates |22| . 
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Appendix A. The matching coefficients 

The values of the matching coefficients a defined in chapter 3 read 
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